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Abstract 

Effect of energetic disorder on charge carrier transport in organic materials has been reex- 
amined. A reliable method for mobihty calculation and subsequent evaluation of relevant 
disorder parameters has been discussed. This method is well suited for a direct calculation 
of the magnitude of dipolar disorder cidip in polar organic materials from the current tran- 
sients. Calculation of (Xdip for several transport materials with varying concentration of polar 
dopants gives concentration dependences that are in reasonable agreement with theoretical 
predictions. A possible solution of the puzzle concerning the disorder effect on the mobility 
temperature dependence has been suggested. 

Keywords: electron transport, drift mobility, energetic disorder 

1 Introduction 

Modern paradigm of the hopping charge transport in disordered organic materials (molecu- 
larly doped polymers, low-molecular weight organic glasses, amorphous polymers) states that 
the most important factor governing the behavior of a charge carrier is energetic disorder 
[H [2] . It is generally accepted that hopping of a charge carrier in random energy landscape 
U{r) describes very well all major features observed in transport experiments. Crucial ingre- 
dients of the most successful realization of this idea includes the Gaussian density of localized 
transport states 
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and Miller- Abraham hopping rate [3] ; here A^o is a total concentration of transport sites and 
a is rms disorder. If we assume no spatial correlation between energies of different trans- 
port sites (a famous Gaussian Disorder Model (GDM) [1]), then the field and temperature 
dependence of the carrier drift mobility ^ could be described as [1] 



fj, = fio exp (— 2aa) exp 
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, a = aP, P=l/kT, (2) 



where a is a lattice scale for the equivalent homogeneous lattice model with the same con- 
centration of transport sites (a scales as c~^/^ with concentration c of transport sites), a is 
an inverse localization radius of the wave function of the transport level, parameter C scales 
as a and parameter S describes the positional disorder (more precisely, if S < 1.5, then the 
term in eq [2] should be replaced by 2.25). 

In fact, the Poole- Frenkel (PF) mobility field dependence In/i oc V^,usually observed 
in experiments, could be described by the GDM only in rather narrow field range around 
E ^ 1 X 10^ V/cm. This limitation was naturally removed in the more advanced model of 
correlated disorder [H El El E]- Indeed, one of the basic assumptions of the GDM, i.e. the 
assumption of the absence of spatial correlations in the random energy landscape proves to 
be spectacularly wrong in organic materials. In these materials interaction of a charge carrier 
with randomly located and oriented dipoles (the model of dipolar glass [1]) or quadrupoles 
(the model of quadrupolar glass [8]) generates highly correlated energy landscape with corre- 
lation function C(r) = {U{f)U{0)) decaying as 1/r for the dipolar glass [HE] and as for 
the quadrupolar glass [8]. The model of dipolar glass (DG) is a good candidate to describe 
transport properties of polar disordered organic materials, while the model of quadrupolar 
glass (QG) naturally describes non-polar materials having zero dipolar moment but nonzero 
quadrupolar moment; quadrupolar moments are sufficiently high for transport molecules 
having local polar groups with compensating dipolar moments [H]. 

Correlated landscape in the dipolar glass naturally produces the Poole-Frenkel mobility 
field dependence, as it was shown in one dimensional (ID) model of charge transport, where 
correlation function with power-law decay C{r) oc 1/r" leads to the mobility field dependence 
In /i oc [6] . This conclusion has been confirmed by the 3D Monte Carlo simulation 

[7j, where the mobility obeys the relation 



/i = /io exp (— 2aa) exp 
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with C ~ 0.78 and P ^ 2 for the case of totally filled lattice. In striking contrast with the 
GDM, the model of dipolar glass (sometimes called the correlated disorder model) demon- 



2 



strates the Poole-Frenkel field dependence in the field range from lO'^ V/cm to 10^ V/cm, 
in good agreement with experimental data. For nonpolar materials the QG model gives 
In/i oc E^^^ |8], but this dependence is hardly distinguishable from the PF dependence 
in rather narrow field range typical for time-of-fiight (TOF) data in nonpolar materials 

Pini [ID- 
Recent paper by Schein and Tyutnev [12] put forward a challenge to this picture. They 
analyzed a vast set of data on the temperature dependence of the mobility, extrapolated to 
the zero electric field (in short, mobility temperature dependence). According to eqs[2]and 
[3l In ii{E ^ 0) oc — 0"^, giving us an obvious opportunity to calculate a from the TOF data. 
If we assume that different sources of energetic disorder are independent, then the total o"^ 
is a sum of contributions from the individual sources 



One of the most important sources of energetic disorder in polar organic materials is a dipolar 
disorder, properly described by the DG model. For a lattice version of the DG model on a 
simple cubic lattice 

here c is a fraction of sites, occupied by dipoles, p is a dipole moment, and e is a dielectric 
constant of the medium [TSj [H]. 

We should expect that other sources provide contributions to the sum in eqlH e.g. charge- 
quadrupolar interactions or van-der-Waals (charge- induced dipole) interactions. Nonethe- 
less, if we consider highly polar transport dopants (with p = 3 — 5 D) in nonpolar polymer 
binder, e.g. polystyrene (PS), then we could reasonably expect that the dipolar term (eq 
|5]) gives a dominant contribution to the total sum, so a should decrease with the decrease 
of the concentration of transport dopant c. Yet in |12j it was noted that experimental data 
demonstrate approximately constant magnitude of a for c varying from 10% to 70%. The 
only way to provide an agreement between eq [5] and experimental observations seems to 
suggest a contribution to eq H] which increases with the decrease of c. Although such con- 
tributions do really exist (for example, charge-quadrupolar or van-der-Waals interactions of 
charge carrier with polystyrene units), their contributions are expected to be much smaller 
than the dominant dipolar contribution. According to the analysis of Schein and Tyutnev, 
the situation is very general. For this reason they concluded that the most possible reason 
for the mobility temperature dependence is a contribution from some intramolecular mech- 
anism. A potential reader of their paper may ask an unavoidable question: where is the 
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disorder? Why does not it manifest itself in the temperature dependence of the mobihty? In 
this paper we are going to discuss this fundamental problem affecting any possible disorder 
model of charge carrier transport in organic materials. 



2 How to calculate (Jdip 



In order to resolve a problem of the mysterious absence of the disorder effect, let us start 
with the examination of the crucial difference between the GDM and correlated models. 
We already noted that the PF behavior cannot be described in any convincing way by the 
GDM but arises naturally in the DC model. The most important transport property of the 
correlated disorder is a relation between the spatial decay of the correlation function and 
the mobility field dependence [6j. This fundamental relation offers a possibility to calculate 
a contribution from the dipolar disorder to the total a directly from the TOP data. 

This can be done because according to the correlated disorder formalism (the simplest 
tool here is a ID transport model developed in ^) the most important contribution to the 
mobility field dependence for moderate field is provided by the most correlated component 
of the disorder, i.e. the component with the slowest decay of the correlation function C(r). 
For example, for the algebraic correlation function C{r) = Aa"^ {O'/'^Y most important 

case of strong disorder cr/3 ^ 1 the ID model gives 

n 

/ 1 \ 1 f eaE\ "+i 

\n^i/^iQ = G{T,E)^-2aa-a'^(3'^+{l + -\a(3{Ana(3)~ . (6) 

As we already noted, n = 1 for the dipolar disorder and n = 3 for the quadrupolar disorder; in 
both cases a dimensionless constant A ~ 1 [UIH]. Formally, the GDM case may be considered 
as the limit n ^ oo. Mobility field dependence can be estimated by the derivative 

- = -(A„.«.«(— ) . (7) 

Taking into account that a (3 ^ 1 and eaE/a < 1, we see that the right hand side in 
eq [7] decreases with the increase of n, i.e. with the decrease of the strength of spatial 
correlation (contribution of the multiplier n^+r is unimportant). Thus, eq [7] indicates that 
the most correlated disorder provides the most important contribution to the mobility field 
dependence. 

Dominant role of the dipolar disorder could be easily demonstrated for the case of the 
two component disorder 

t/(rO = f/dip(f) + f/„(f), (8) 
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where some less correlated independent random energy Un{r) gives an additional contribution 
to the total disorder. If the correlation function for this contribution has a form C„(r) = 
AnO'n^ with n > 1, then the ID model shows that at the low field boundary of the PF 
region eaf3E > ^ the additional contribution to the mobihty field dependence is 

small, if 



0": 



n ■ 



-«1 



(9) 



(see Appendix). In many organic materials PcTdip ^ 1 and dipolar contribution dominates 
the mobility field dependence even if o"„ ~ o"dip. 

The result of the ID model could be supported by the direct 3D Monte Carlo simulation 
of charge carrier transport in the energy landscape comprised by the mixture of the dipolar 
and non-correlated GDM disorder (see Figure [1]). In Figure the decrease of the slope 
S = 0^/2 of the mobility field dependence (or decrease of |^ in Figure [T)d) for eaE/a < 1 
with the decrease of dipolar contribution (Xdip to the total a is clearly visible. 
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Figure 1: Mobility field dependences for a mixture of dipolar and non-correlated GDM 
disorder for the constant value of the total a; here [ig = A^oexp (— 2aa), a"^ = a^jp + ctq^j^, 
and the ratio r = c^jp/cr^ is indicated at the corresponding curve. Transients have been 
simulated for kT/a = 0.26. U a ^ 1 nm and a ^ 0.1 eV, then eaE/a ^ 1 ioi E = 1 x 10^ 
V/cm. Note the transformation of the mobility field dependence for moderate field from the 
PF type for r = 1 to the linear dependence InficcE for r = 0. This very dependence is an 
intrinsic property of the GDM [7]. 
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Hence, the temperature dependence of the slope S(T) 



siT) = ^ = ^ i-'^' - r) (10) 

in polar materials could be used for a direct evaluation of (Xdip- Other sources of disorder, 
as well as the contribution from traps [15], provide much less significant contributions to 
S{T) for eaE < cTdip- Situation is especially favorable for highly polar dopants in PS matrix 
because the dipolar moment of styrene molecule is only 0.4 D [16j, and we can expect that the 
quadrupolar moment of styrene is rather low too (styrene molecule is built only by carbon 
and hydrogen atoms and does not contain polar groups). Traps do not affect S(T), either 
[T7] . In future we will use notations as for the rms disorder, calculated from S{T), and ax for 
the corresponding rms disorder, calculated from the temperature dependence of fi{E 0) . 
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Figure 2: Dependence of a on the concentration of polar constituents [T8l IT6]: as (empty 
symbols) and ax (filled symbols), correspondingly. Squares correspond to the system 
DEH:PS [IH], circles correspond to the system with copolymer PS-BA as a binder [16], 
and lines are shown as a guide for an eye. Data from the paper |l9j are not shown to avoid 
overcrowding of the figure. 

Information we needed is pretty rare in literature and we found only four papers where 
the data for S{T) have been provided for different concentrations of polar constituents 
[HI [HI [20l [19]. Fortunately, these papers describe very different situations. Two papers 
consider dependence of the mobility on the concentration of highly polar transport dopants 
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DEH (4-diethylaminobenzaldehyde diphenylhydrazon, p = 3.61 D) [T8] and TPM-E (triph- 
enylmethane derivative, p = 2.1 D) [19j in PS matrix. The third paper describes dependence 
of the mobihty on the concentration of polar units in the copolymer of styrene with polar 
monomer BA (butylacrilate, p = 1.52 D) [16], serving as a polymer binder, and the last 
paper considers an effect of the inert polar dopant TAP (t-amylphthalonitrile, p = 6.6 D) 
[20] on the charge transport in weakly polar transport material. In the two last cases concen- 
tration of the transport dopant was fixed. In all cases calculated cxdip demonstrates the same 
expected trend (see Figures [2] and [3]) : it decreases with the decrease of the concentration 
of polar constituents with the expected magnitude of the decrease (apart from the case of 
TPM-E:PS, where the decrease is rather small). 
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Figure 3: Dependence of cr^ on the concentration of inert polar dopant TAP [20]: cr| (empty 
squares) and erf. (filled squares), correspondingly. Lines are shown as a guide for an eye. 

It is worth noting the almost linear dependence of (t| on c for the transport material 
with small concentration of polar dopant TAP, in good agreement with eq[5] (Figure [3]). We 
should expect that the model of ideal dipolar glass without any orientational correlation 
between dipoles is best suited for a description of this very situation: low concentration of 
polar dopant in amorphous material. Concentration dependencies of as in Figure [2] do not 
obey eq[5l This deviation may reflect a partial orientational ordering of dipolar molecules for 
high concentration of dipoles or a defect of the particular method of the mobility calculation, 
discussed in detail in the next section. 
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We would like to emphasize that the very scarcity of the published data on the tem- 
perature dependence of 5* could be considered as a strong argument in favor of general 
agreement between predictions of the correlated disorder model and experiment. Four pa- 
pers p!6l ITSl l20l IT9] may be regarded as four random picks chosen from the vast variety of 
the TOF data. Remarkably, these four random picks immediately demonstrate reasonable 
agreement with theoretical predictions. 

If the dipole moments of all constituents of the transport layer are low, then the DG 
model is clearly inadequate for the description of the charge transport. The most promising 
candidate in that case is the QG model. Here we provide only a very brief summary of the 
Monte Carlo simulation results for the QG model, more complete report will be published 
elsewhere. General approach is almost identical to the one described in [7j. Results of the 
simulation (see Figure H]) could be summarized as 



fi = /ioexp (— 2aa) exp —0.37a^ + Cq (a^^^ — Tq) {eaE/a 



3/4 



fir 



with Cq = 0.87 and Tq = 1.91. In nonpolar organic materials mobility field curves should 
be fitted to eq [HI and the corresponding quadrupolar cr should be estimated from the 
temperature dependence of the slope of the dependence In n vs 



3 Dangers of low field mobility 

In the previous section it was demonstrated that dipolar energetic disorder in organic mate- 
rials can be directly evaluated from the TOF data. Result of this evaluation demonstrates 
a reasonable behavior of a^ip with respect to concentration of polar constituents. Magni- 
tude of (Jdip (and its variation with concentration) again is in reasonable agreement with the 
theoretical estimation by eqO Now we can reformulate the crucial question raised in the In- 
troduction in a different way: why does the energetic disorder manifest itself in the mobility 
field dependence, but not in the mobility temperature dependence? Why the difference? 

Discussion in the previous section unambiguously demonstrates that disorder does affect 
hopping transport of charge carriers in organic materials. We cannot consider as a possible 
alternative explanation an assumption that some totally unknown mechanism is responsible 
for the PF field dependence and, at the same time, it just by accident provides the same 
dependence of transport parameters on the concentration of polar dopants as a well known 
mechanism of the correlated dipolar disorder. 

If disorder does affect charge transport, then there is no possibility to avoid a corre- 
sponding term proportional to 0"^; in the mobility temperature dependence. According to 
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Figure 4: Mobility field dependence in the QG model for different values of kT/a (from the 
top curve downward); straight lines indicate the fit for eq[TTl Plot of the simulation data 
in the usual PF presentation In jj, vs E^^"^ demonstrates much stronger deviation from the 
linearity in the weak field region. 

the estimations of eqs O and [3], this term is not negligible but quite comparable to the total 
cr^. Hence, only two possibilities are feasible: either some compensating mechanism provides 
a contribution to the total a (more rigorously, to the slope of the mobility temperature de- 
pendence) or some defects of the experimental procedure masquerade the true dependence 
of a on c. Let us consider both possibilities. 



3.1 Possible physical compensating mechanisms 

In the previous sections we demonstrated that S{T) gives a robust estimation for the mag- 
nitude of (Tdip. Temperature dependence of the low field mobility is not so easily tractable. 
A host of various contributions adds to the resultant dependence: charge-dipole, charge- 
quadrupole, van der Waals (charge-induced dipole) interactions, polaron contribution of a 
simple Arrhenius type t21]|, and the similar contribution from traps. Reliable separation of 
the individual contributions to the resultant temperature dependence is impossible. As it 
was already mentioned, some contributions could provide compensating effects, for exam- 
ple contributions originated from the carrier interactions with polymer binder, though the 
magnitude of the effect is expected to be insufficient for a total compensation in the case 
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of highly polar transport dopant in the PS binder. There is another possibility for compen- 
sation because for a low concentration of the dopant charge transport should become more 
sensitive to the energetic disorder: exponential decrease of the hopping rate for large average 
distances between dopant molecules results in carrier hopping to the nearest neighbor with- 
out regards to its energy. This means that the effective ut should increase with the decrease 
of c, but this effect is hardly very important for c > 10 — 15%. 

Effect of impurities (possible traps) deserves a special discussion. First of all, traps do 
affect effective magnitude of ut, at least the one estimated using the GDM [22]. Effect of 
impurities could be responsible for different values of (Tt observed for the same material. 
For example, dependence of cr-r on c for the same system DEH:PS was found to be different 
in different papers [HI |23]. In the earlier paper [23] a constant value ctt ~ 0.13 eV was 
found. In the later paper [18] more rigorous purifying procedures were used and very careful 
preparation procedure was undertaken. As a result, variation ^ 15% in the magnitude of 
(Tt over studied concentration range has been found (see the corresponding curve in Figure 
|2]). This is certainly beyond the experimental error range. For the same concentration of 
DEH, magnitude of 0"^, reported in Refs. [18] and [23], differs up to 20-25%. This means 
that the slope of the mobility temperature dependence could differ up to 40% for different 
procedures of material purification and sample preparation. In [12] the result of [18] is 
mistakenly regarded as corroborating the conception of intramolecular origin of 0"^, while in 
fact comparison of the results of Refs. |18j and [23] suggests that more careful examination 
of the trap effect on ctt, calculated from the TOF data, is needed. 

Some details of the preparation procedure, exploited in [18], are unusual. For example, 
coating of the transport layer was performed under the red light to avoid photochemical 
creation of traps. Quite probably, this very careful preparation procedure gave a possibility 
to avoid excessive amount of traps and obtain more reliable values of ax- 

Note, that polymers are difficult to purify; to some extent, all of them contain vari- 
ous impurities, possibly serving as traps. Analysis of charge transport in trap-containing 
materials shows that increase in trap concentration could be interpreted as an increase of 
the effective cr-r [24j. In molecularly doped polymers decrease of the dopant concentration 
means a simultaneous increase of the polymer concentration and, thus, possible increase in 
concentration of traps. This could provide an additional compensating mechanism. 

3.2 Deficiency of the mobility calculation 

So far we considered a possible compensation that is relayted to some intrinsic physical 
mechanism. Yet, quite probably, the most important reason for the compensation is a 
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peculiar procedure for the calculation of the drift mobility, employed in the majority of 
experimental papers. 

There are two most frequently used procedures for the calculation of fi = L/Et from 
the TOF data in double linear current vs time presentation: 1) calculation that uses time ti 
determined by the intersection of asymptotes to the plateau and trailing edge of the transient 
and 2) calculation that uses time ti/2 for current to reach half of its plateau value; here L is 
a thickness of the transport layer. The first method is a method of choice for most papers. 
One can find in literature the statement that the difference between two procedures is not 
very important for determination of the temperature and field dependence of the mobility 
[251 EH]. This is not true. 
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Figure 5: Mobility field dependence for two methods of the mobility calculation: /ij (filled 
squares) and /ii/2 (empty squares) [27]. Straight lines indicate best fits. 

Calculation, based on ti, significantly overestimates contribution of fast carriers. Defi- 
ciency of this calculation is especially evident in the case of materials with traps. In trap-free 
materials the use of /Xj leads only to the change of the slope of the mobility field dependence. 
In the presence of traps the use of fii sometimes distorts even the functional type of the 
mobility field dependence, bringing it to the linear one In/Xj oc E. This behavior was ob- 
served in experiment [28], but it is an artifact, directly related to the inferior method of the 
mobility calculation [29j. Procedure that uses ti/2 is much more reliable and generally agrees 
well with the standard definition of the mobility as n = v/E, where v is an average carrier 
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velocity. 

Unfortunately, errors, related to the use of tj, are not limited to the mobility field de- 
pendence, they also affect the mobility temperature dependence. As an illustration we pro- 
vide Figure [5], where the temperature dependence of fii and yUi/2 for tri-]9-tolylamine doped 
polystyrene is shown. Slopes of the dependencies differ by approximately 20%, thus the 
values of ctt, calculated from the slopes, do differ too. Similar difference was found for the 
data presented in [30] • Again, as in the case of the temperature dependence of S, the data 
for raw transients (note that here we need transients for several values of T) is very difficult 
to find in literature. 

It is very important to understand the full significance of the described fact. If we cal- 
culate (Tt from the slope of the mobility temperature dependence for different values of c 
and use the mobility /ij, then there is absolutely no guarantee that its relative difference 
from the value of ax, calculated from fii/2, remains constant in the whole range of concen- 
tration. Remember, /ii/2 is a much better approximation for a true mobility. Hence, the 
obtained concentration dependence of a-r is certainly not very reliable and, quite probably, 
has considerable errors. 

3.3 Non-dispersive or dispersive transport? 

There is another aspect of the mobility calculation procedure that can be a possible source 
of errors. It was claimed in [12] that for all cases used in the analysis of charge transport, 
transport regimes are essentially non-dispersive. This is not true. For example, Borsenberger 
et al. explicitly emphasized that transients for highly polar dopant DTNA (di-;?-tolyl-p- 
nitrophenilamine, p = 5.78 D) in PS are highly dispersive and no plateau was observed in 
the whole field, temperature and concentration range; sometimes the transit time has been 
even determined by the intersection of the asymptotes in double logarithmic current vs time 
representation [31]. For this very reason ax was even calculated by the formula 

ln/i//io ~ -2aa - (12) 

that differs from eqO It was suggested by Bassler and Borsenberger [32] for the treatment 
of low field mobility in the case of dispersive transport. Particular material DTNA:PS could 
be an exceptional case in the sense that in the whole range of all relevant parameters c, T, 
and E transients demonstrate no visible plateau, but for low temperature (T < 230 K) this 
kind of transients is generally a rule and not an exception [321 [331 [3l] ; how the transit time ti 
(or ^1/2) could be reliably calculated in double linear plot is not very clear in this situation. 
In fact, sometimes the problem of the choice between non- dispersive and dispersive trans- 
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port regime is even more tricky than described in the previous paragraph. For example, let 
us consider again the case of DTNA-doped PS. In the original paper [31] transients were con- 
sidered as dispersive because no visible plateau was observed. Yet it is very well known that 
non-dispersive transport (described by the usual diffusion equation with well defined v and 
diffusion coefficient D) can routinely produce the same kind of transients if vL/ D < 5 — 10. 
Such transients have no well defined plateau but could be easily discriminated from the true 
dispersive transients by their asymptotes for t — oo. True dispersive transport transients 
obey the law I{t) oc |35], while for a normal diffusion I{t) oc exp {—v'^t/AD), putting 

aside an unimportant power-law factor. 

This test immediately demonstrates that at least some transients, provided in [31], are 
not dispersive: the transient can be fitted with reasonable accuracy by the solution of the 
diffusion equation [36l [371 EH] 



/(t) oc 




vH\ f [L-vtf 

exp I — — exp ' 



1 +v erf ^ ] +erf { -^L=\ 



AD J ' \ ADt 

(13) 

where erf(a;) is the error function, and the large time asymptotic is of the exponential kind 
(see Figure [6j). We can safely conclude that in some cases charge transport in DTNA:PS 
occurs in the non- dispersive regime, and yet it was always analyzed using eq[T21 suitable only 
for highly dispersive regime! All this confusion means that the true magnitude of ax and its 
dependence on c are very different from provided in [31]. One can learn a very important 
lesson from this example: in order to calculate mobility from transients without plateau the 
very first step should be an attempt to fit the transient to eq[T3l 

One can conclude that the data, discussed in [12], had been obtained both for non- 
dispersive and dispersive regimes. Even worse, in some cases a non- dispersive mobility had 
been treated as a dispersive one, leading to the patently wrong calculation of ax- We believe 
that in this situation it is impossible to reveal the true dependence of the rms disorder on 
the concentration of polar constituents. 



3.4 Possible hints for improper treatment of TOF data 

Discussion of possible difficulties relevant to the analysis of TOF transients and evaluation of 
the disorder parameters from the TOF data should be finished with one additional remark. 
Although current transients are of primary importance for a detection of possible deficiencies 
of the mobility calculation procedure, sometimes the hints for the improper treatment of the 
transport data could be obtained directly from fi{E) or /i(T) curves. 

One very well known example is an artifact of the mobility decrease with the increase of 
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Figure 6: Transport transient from Figure 2 of [31J (crosses) and the best fit for eq[T^ (line). 
Note that ln/(t) oc t for t > L/v. 

E in weak field region. Earlier this phenomenon was attributed to the effect of field induced 
traps [T] , though this explanation for weak field region always looks pretty suspicious. Later 
it was found that the true reason for this strange behavior is a deficiency in the mobility 
calculation procedure, i.e. the use of /Xj [361 EH [39] • If fitting of the transients to eg [T3] 
is used, then the mobility curve becomes monotonously increasing with E. The reason 
for the unexpected increase of /i in weak field region is a trivial contribution from the 
normal diffusion, described by the diffusion coefficient D. It is worth noting that the fitting 
procedure, suggested by Hirao et al. [MIET], produces an additional important result: it 
removes the spurious dependence of ii on the thickness of the transport layer L {D was also 
found to be independent of L). 

Artifact of the mobility field dependence in weak field is directly related to the main 
subject of our consideration. Mobility that decays with the increase of E was reported for 
transport material TAPC:PC (l,l-bis-(di-4-tolylaminophenyl)cyclohexane in polycarbonate) 
doped with small amounts of inert polar dopants o-, m-, and p-DNB (isomers of dinitroben- 
zene) [ID]. For m-DNB and p-DNB the mobility decays with the increase of E in the whole 
field range (up to ~ 1 x 10^ V/cm). Such behavior for room (or lower) temperature is 
extremely unusual and, very probably, indicates an improper method of the mobility calcula- 
tion from the TOF data. Correspondingly, magnitude of ctt, provided in [lO] and extensively 
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discussed by Schein and Tyutnev [12], is hardly very reliable. 

Highly unusual value of some relevant parameter, obtained from the interpretation of 
TOF data, could serve as another alarm signal. As an example of such alarm signal we 
consider an unusual value of the radius of the localization Rq of the wave function of transport 
site, obtained for DEH:PS [23j. It was found that Rq = 1/a = 4.8 A, while for all other 
transport materials Rq typically falls in the range 1 — 2.5 A [HI HJl |l3l HI]. We believe that 
this is another indication that something goes wrong with the mobility calculation; again, 
the use of /ij is a probable culprit. 

4 Checklist for the reliable mobility calculation and 
subsequent evaluation of the disorder parameters 

In this section we would like to emphasize principal details of the more safe and reliable 
procedure for the analysis of the TOF data in order to elucidate disorder parameters. 

• First of all, wherever possible, transport should be analyzed for a non- dispersive regime 
only. In this regime transport parameters are more sensitive for peculiarities of the 
random energy landscape (comparison of eqs [2] and [12] supports this statement for the 
GDM). 

• If a well defined plateau of the transient is observed, then the drift mobility may be 
calculated using ti/2. Use of ti should be completely eliminated from the experimental 
practice. 

• If there is no well defined plateau, then the transient should be fitted to eg [T3l according 
to the procedure, suggested by Hirao et al. [36j. Only in the case, when the quality 
of the fit is clearly inadequate, transport should be analyzed in the framework of the 
dispersive regime. 

• For the case of dipolar disorder, adip should be calculated from the temperature de- 
pendence of the slope of the mobility field dependence S{T) = for moderate field 
using eqO The use of the DG model instead of the GDM is important, because the 
very essence of this approach is based on the transport properties of correlated en- 
ergy landscape. This method is almost completely insensitive to the contribution from 
traps, polaron effects, and van der Waals interaction. The quadrupolar contribution 
probably could not be totally eliminated in this way, but we may reasonably expect 
that it is not dominant for polar materials. 
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• For nonpolar materials the magnitude of quadrupolar disorder should be calculated 
from the temperature dependence of the corresponding slope Sq{T) = ^p/j according 
to eq[IIl Unfortunately, description of the mobility field and temperature dependence 
for real organic materials with contributions from different sources of disorder having 
different correlation properties cannot be described by any simple unified formula. 

Discussion of the dependence of ctt on the concentration of dipolar dopants and other 
relevant parameters should be performed only after calculation of fi and as in the proper 
way. We hope that the realization of this program will provide us with much more reliable 
knowledge of the crucial parameters of the energy landscape in disordered organic materials. 

5 Conclusion 

This paper is a direct response to the recent review by Schein and Tyutnev fi2\. They 
analyzed a vast set of available transport data for disordered organic materials and came 
to the conclusion that in some mysterious way energetic disorder is not important for the 
mobility temperature dependence. The strongest argument in support of this point of view 
is the independence of disorder magnitude a on the concentration of highly polar transport 
dopants. This is very unexpected because the dipolar contribution to the disorder should 
depend on the concentration of polar dopants. Schein and Tyutnev suggested that some 
intramolecular mechanism is responsible for the mobility temperature dependence. 

We have analyzed essentially the same set of data but in a different manner and came 
to a very different conclusion. Our analysis shows that the dipolar contribution could be 
directly calculated from the temperature dependence of the slope S of the mobility field 
dependence. Results of such calculations immediately demonstrate a reasonable agreement 
between experiment and theory, though available experimental data is extremely limited. 

We have shown that a very popular method of the mobility calculation by the time of 
intersection of asymptotes to the plateau region and trailing edge of the transient is a possible 
reason why the calculations of a give misleading results. A simple yet accurate procedure 
of the mobility calculation and subsequent analysis of the mobility field and temperature 
dependences for evaluation of the disorder parameters has been discussed. 

We would like to emphasize that a very significant contribution to the realization of the 
suggested program could be fulfilled by a simple recalculation of the mobility from already 
obtained transients (if available) and subsequent analysis of the mobility field dependence 
according to the lines suggested in the previous section. Also we would like to note that 
a direct access to the raw TOF data (current transients) should significantly improve our 
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ability to check consistency between experimental data and various theoretical models and 
would be considered as God's gift by all theoreticians. Data for TOF transients could be 
provided as supplements to the experimental papers, this possibihty is already provided by 
many leading scientific journals. 
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Appendix. Mixed disorder in ID model 

Suppose that the total random energy U (r) is a sum of two independent terms 

t/(r) = f/dip(r) + f/n(f), (14) 



where Udip{r) is a dipolar contribution with the correlation function Cdip(r) = ^dipO"dip^ 
while the second term provides an additional energetic disorder with the correlation function 
Cn(^ = ^nO'n^, and n > 1. For independent disorder contributions the resulting correlation 
function is the sum of two terms 

C(f) = Cdip(f) + C„(f). (15) 

In the ID model the mobility is 

Ej^ dxexp {-e(3Ex + [0(0] - C(a;)]} ' ^^^^ 

where Dq is a bare (microscopic) diffusion coefficient [HlllS]. For the case of strong disorder 
l3a ^ 1 a suitable tool for the calculation of integral [16] is a saddle point approximation. For 
example, eq[6]is a result of this approximation. In fact, for a pure dipolar disorder the PF 
dependence arises exactly in the case when this approximation is valid. In the saddle point 
approximation 

/i oc exp [-C(0)/?2 + R{xs)] , R{x) = ex(3E + C(x)/3^ C{0) = + a^, (17) 
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and solution of the equation 



^ = e0E- A,,,fia%^^^ - nA^pV„^ = 0. (18) 

Let us find when the mobihty field dependence retains the PF form. This is the case when 
the third term in eq [18] is a small correction. If we retain the first two terms in eq [181 then 
the solution is 

-s=aZ = a^--^) , (19) 



and the third term is a correction if 



nAn al 1 



Saddle point approximation is valid if three conditions are valid 



.3 d'R 



< 1, 6xs < x° > a, (21) 



where 5x,^ = ( ^ 



-1/2 



. The first condition means that the Gaussian approximation 

is valid in the vicinity of the minimum of R{x), the second one means that the vicinity of 
the minimum does not reach the boundary of the integration domain, and the third one 
means that we can safely use power- law asymptotic for the correlation functions C{x). All 
inequalities could be condensed as 

1 « Z « (3^al^, (22) 



and here and in future we drop all coefficients such as 3/-\/2, A^ip, or An. This double 
inequahty is equivalent to [P'^o'dip) ^ ^ ea^E <^ P^o'dip- The low field boundary of the 
PF region can be estimated as Zi < P'^cr^ip. Substituting Z with Zi in eq [20] we obtain 
eq [9] as a final necessary condition for the small contribution of the additional non-dipolar 
disorder in the low PF region ea(3E > (o"^ip/3^) ^. The same estimation could be obtained by 
comparison of the dominant contribution to the mobility field dependence with the additional 
term oc AnCr^- 
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